The ion-induced charge- exchange X-ray emission of the Jovian 
Auroras: Magnetospheric or solar wind origin? 

Yawei Hui and David R. Schultz 

Q\- Physics Division, Oak Ridge National Laboratory 

o : 

Bldg. 6010, Oak Ridge, TN 37831 USA 
huiyOornl.gov, schultzd@ornl.gov 
^ . Vasili A. Kharchenko 

■ Physics Department, University of Connecticut 

■ 2152 Hillside Road, U-3046, Storrs, CT 06269 USA 
q' kharchenko@phys.uconn.edu 

a: 

Phillip C. Standi 

cd . 

Department of Physics and Astronomy and the Center for Simulational Physics, University 

■ of Georgia 

<n : 

!> ' Athens, GA 30602 USA 

stancil@physast.uga.edu 

o : 

Q\ . Thomas E. Cravens 

O ■ 

J> ■ Department of Physics and Astronomy, University of Kansas 

^ ■ 1251 Wescoe Hall Drive, Lawrence, KS 66045 USA 

cravensOku . edu 
Carey M. Lisse 
Johns Hopkins University Applied Physics Laboratory 
SD/SRE, MP3-E167, 11100 Johns Hopkins Road, Laurel MD 20723 USA 

carey . lisse@jhuapl . edu 
and 



-2- 



Alexander Dalgarno 
Harvard- Smithsonian Center for Astrophysics 
60 Garden Street, Cambridge, Massachusetts 02138 USA 
adalgarnoOcf a . harvard . edu 

ABSTRACT 

A new and more comprehensive model of charge-exchange induced X-ray 
emission, due to ions precipitating into the Jovian atmosphere near the poles, 
has been used to analyze spectral observations made by the Chandra X-ray Ob- 
servatory. The model includes for the first time carbon ions, in addition to the 
oxygen and sulfur ions previously considered, in order to account for possible ion 
origins from both the solar wind and the Jovian magnetosphere. By comparing 
the model spectra with newly reprocessed Chandra observations, we conclude 
that carbon ion emission provides a negligible contribution, suggesting that solar 
wind ions are not responsible for the observed polar X-rays. In addition, results 
of the model fits to observations support the previously estimated seeding kinetic 
energies of the precipitating ions (~0.7-2 MeV/u), but infer a different relative 
sulfur to oxygen abundance ratio for these Chandra observations. 

Subject headings: atomic processes — planets and satellites: individual (Jupiter) 
- X-rays: individual (Jupiter) 



1. Introduction 



Since the first generation of X-ray observatories, Jupiter has been one of the primary 
solar system objects of interest. In fact, Jovian X-ray emission was fir st detected by the 



Einstein observatory (IMetzger et al.lll9831 ). then studied with ROSAT (IWaite et al.l Il994 



Gladstone et al.lll998l ). and two distinct components were identified, namely emission from 
the high-latitude (polar) and low-latitude (equatorial) regions. Various explanations of the 
origin of these emissions were put forth but could not be stringently tested due to the poor 
spatial and spectral (S/AE < 1) resolution of the early measurements. However, recent ob- 
servatio ns with Chandra and XMM-Newton have provided unprecedented resolution on t he 
planet Jdadstone et al.ll2002l : felsner et al.lbood : iBranduardi-Ravmont et alibooi . l2007al lbl. 
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20081 ; iBhardwaj et al.ll2005l . 120061 ) and thorough analysis of the observational data have con- 
firmed and refined the characteristic difference between Jovian polar and disk emission. It 
is now generally accepted that the disk component is due to the scattering and fluores- 



cence of solar X-rays in the atmosphere (IMaurellis et al.l I200CH ; iBranduardi-Raymont et al. 



2007bl : iBhardwai et aDbood . bood : bravens et al-lbood ). In contrast, a mounting body of 



evidence favors the charge-exchange (CX) model as the basic explanation of the Jovian 
X-ray polar auroras. This model assumes highly energetic heavy ions precipitating into 



cited states ( 


Metzger et al. 


1983; 


Horanvi et al. 


1988:IWaite et al. 


1994 


; Cravens et al. 


1995 


Kharchenko et al. 


1998|; 


Liu & Schultz 


1999; 


Kharchenko et al. 


2006. 2 


008 


)■ 



However, fundamental questions remain unanswered, such as What is the detailed mech- 
anism of ion acceleration before emission; What is the source of the heavy ions (do they 
originate from the Jovian magnetosphere and are therefore closely coupled with its ion- 
rich satellites, or are they of solar wind origin); and What are the composi tion of the ion 



fluxe s and their initial kinetic energies at the top of the Jovian atmosphere (ICravens et al. 



20031 1? Observations and interpretive models answering these questions would advance our 
knowledge of the Jovian system and, by extension, of other planetary atmospheres and en- 
vironments. Toward these ends, many studies have been carried out seeking to explain the 
observed spectral features by fitting them with a continu um plus ion emission lines (e.g., 
Eisner et all 1200,4 iBranduardi-Ravmont et al.lbooi l2007al lbl). In order to reflect the intrin- 
sic correlations between the intensities of different X-ray lines induced by individual ions, 
CX models were developed to produce synthetic spectra. These spectra could, in turn, be 
used to infer parameters such as relative ion abundances and initial energies required to 
account for the obse r vations. Early CX models were restricted to oxygen ion precipitation 
( 1 Cravens et al.l Il995l ; iKharchenko et al.l Il998l ; iLiu fc Schultz] 119991 ) owing to the lack of ex- 
tensive atomic data required by the models (i.e., energy- dependent, state-selective charge 
exchange cross sections for all ionization stages of the ions and transition probabilities for 
radiative and non-radiative de-excitation, etc.). 

Since oxygen ions contributing to the Jovian X-ray auroras could originate either in the 
vicinity of Jupiter or from the solar wind, the next level of discrimination between these two 
possible sources was provided by including sulfur ions because of their great magnetospheric 
abundance. This was enabled by the advent of larg e scale atomic collision and s tructure 
calculations, and recently included in the CX model (IKharchenko et al.ll2006l . 120081 ). These 
results provided an improved fit to the observations and gave indications of the range of initial 
energies of the ions at the top of Jupiter's atmosphere. Here, we seek to further constrain 
the origin and characteristics of the X-ray emission through (i) the addition of carbon to 
the model, again enabled by newly available atomic data, (ii) significantly increasing the 
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level of detail of the treatment of radiative transitions, and (iii) using a more sophisticated 
methodology to fit the synthetic spectra to observations through systematic variation of the 
input ion abundances and kinetic energies. 

In this Letter, we report the results of this more complete CX model of the Jovian 
auroral X-ray emission. The basic structure and major improvements of the model are 
described in Section [2j In Section [3] a description is given of how our results have b een fit to 
Chandra observations with the X-ray spectral fitting software Xspec (jArnaudlll996l ). Finally, 
discussion and interpretation of the X-ray emission modeling results are given in Section HI 



2. Ion-induced charge-exchange emission model 



A des cription of the bas i c CX model along w i th its development has been given pre 
viou slv bvlKharchenkq et al.l_rtl998h . lLiu fc Schult] Jl999h. IKharchenko fc Dakarnol (boooh. 
and 



2000 



Kharchenko et a 



Bhardwai et al 



(120061. 120081 ). and in various reviews (e.g., iBhardwaj fc Gladstone 
2007^ 

I n brie f, we begin with the Monte Carlo (MC) simulation de- 



veloped by IKharchenko et al.l (120081 ) to treat the ion precipitation and CX collisions during 
the ion's deceleration through the atmosphere. Modeling the energy loss and capture into 
the excited states that seed subsequent radiative emission involves three collision channels 
- charge exchange, target ionization, and electron stripping from the ions. Each channel 
has an associated energy loss contributing to the stopping process, and stripping and charge 
exchange change the ion's charge state. CX leads to population of excited states that result 
in photon emission or, less frequently, decay non-radiatively. A synthetic spectrum is built 
up by tracing all possible cascading paths after these excited states are populated. 

Specifically, we set the total number of ions, N to t, to 1000 in the MC simulation for each 
species considered; increasing this number mak e s neg ligible changes to the final synthetic 
spectrum. As pointed out by IKharchenko et al.l (120081 ). the spectra are not sensitive to the 
initial charge states of the ions for a fixed ion energy at the top of Jupiter's atmosphere, 
so we choose singly charged ions (e.g., + ) to begin the precipitation pr ocess. The energy 



l oss of the ions due to collisions in all three channels are as formulated by IKharchenko et al. 



( 119981 ) . A simulation terminates when the ion slows down to the point that further X-ray 
emission is negligible. 

The most significant improvement in the present work is the incorporation of car- 
bon ions as a new source of X-ray emission, in addition to the previously considered oxy- 
gen and sulfur, in order to introduce an abundant, more distinctly solar wind originating 
species. This is motivated by Chandra and XMM-Newton observations ([Gladstone et al. 
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20021 ; iBranduardi-Raymont et al.ll2008l ) which appear to overturn the early hypothesis that 
the auroral X-ray emission was linked to ions flowing along magnetic field lines from the Io 
plasma torus. These new observations indicated that the X-ray emission comes from lat- 
itudes higher than those which are magnetically connected to flux from the orbital radius 
of Io. Thus, other sources have been su ggested including from more distant r egions of the 
magnetosphere and the solar wind (e.g., iMauk et al.ll2003l ; ICravens et al.l 120031 ) . 



The MC simulation tracks the energy loss and charge state change at each collision 



accor ding to the probabilities established by the relevant atomic cross sections (ISchultz et al. 



20091 ) and, for each CX collision, records the charge state and ion energy. At the end 
of the simulation the CX number distribution, N (q, E), is tabulated, where q and E 
are the charge st ate and the energy , respectively. The calculated probabilities for state- 
selective capture (ISchultz et al.ll2009i ) are then multiplied by N cx (q, E) to obtain the initial 
population X(q,E;j) where 7 represents the energy level (n,£) of each ion populated by a 
CX collision. 

The next step in producing a synthetic spectrum is to calculate the radiative transition 
matrices, T(g;7', 7), which describes the pat h way f rom any initial state (7) to a n y fina l 
state (7'), as described by iKharchenko et al.l (119981 ) and iKharchenko &: Dalgarnol (l2000l ). 
For the present work, w e collected atomic transi tion data from t he Atomic Spectra Database 
( Ralchenko et al. 2008), the Atomic Line List (jvan Hool 19991). and rece nt results for sev- 



eral ion species (IJohnson et al.l |2002| ; Kingston et al 



2002 



Naharll2002f ). In total, we in- 



clude 960 carbon, 954 oxygen, and 1512 sulfur emission lines when calculating the tran- 
sition matrices (but in the subsequent model fit to the X-ray spectra we only consider 
lines with ph oton energy above 200 eV). This is an elaboration of the model beyond our 
earlier work (IKharchenko et al.l 120081 ) which used energy-independent photon yield tables 
and chose only the most prominent lines of oxygen (12 lines) and sulfur (27 lines). Given 
the transition matrices, the photon yields from the cascade process are determined by 
Y(q, E;j', 7) = T(q; r y',j) x X{q,E;^) and the synthetic spectrum is found by summing 
over all ions and charge states. Figure [1] shows an example of a synthetic spectra for the 
case of equal C, O, and S elemental abundances, initial ion charges equal to one, and initial 
ion energies equal to 2 MeV/u. 



3. Spectral Fitting 

With this CX model, we can seek a fit to the observations by varying its inputs as 
in previous work and judging the quality of the resulting comparison. However, in the 
present case, with three chemical species and their differing initial energies, significantly 
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more combinations of the inputs need to be evaluated and differences in the quality of fits 
can be made quantitative rather than qualitative. For example, with the restriction to 
oxygen ions in the earliest models, given the best available atomic collision and transition 
data, the only input to the model was the initial oxygen ion energy. When sulfur atomic data 
was first computed and compiled the inputs were the initial oxygen and sulfur ion energies 
and their relative abundance. Here, we construct an Xspec model that has six independent 
parameters. These are the trial values of the initial energies of the C, O, and S ions at the top 
of the Jovian atmosphere, the relative abundances of C and S with respect to that of O (Ac, 
As), and the emission line width (the full width at half maximum, FWHM, a free parameter 
of the fitting procedure, possibly varying from the observational line width). Because we 
calculate the relative intensities of the emission lines in the CX model, the synthetic spectra 
intensities are given in arbitrary units and the normalization ("NORM" in Tables 1 and 2 
below) required by Xspec absorbs all other factors affecting the flux. 



We have applied this analysis to the Chandra observations made by lElsner et al.l (120051 ) . 
We first reprocessed the original ACIS-S3 data (Observation ID's 3726 and 4418) using 
the newest calibrations and performed all necessary temporal and spatial filtering. New 
algorithms were applied to adjust the pulse height amplitude values in the original Level 1 
event files. The spectra were processed so that each spectral energy bin (wavelength bin) 
had at least ten counts. All the photons were then reprojected into a reference frame fixed 
on the Jovian disk and the spectra were extracted using the latest Chandra data analysis 
package, ciao 4.0.1. 

We attempted fits of the north and south polar observations 3726 and 4418 with C, 
O, and S but found that the carbon ion energy of the fits approached the very low value 
of 0.03 MeV/u. At such energies, the population of highly charged ions producing X-rays 
is negligible. Closely related to this behavior the fits yielded no bound for the carbon to 
oxygen abundance ratios, meaning that the relative amount of carbon required exceeded the 
preset limit we imposed, in this case a C to O ratio of 1000, or dropped to the lower limit, 
zero. To illustrate this behavior of the fits we display in Table [1] results including C and 
omitting it for spectra from the north polar region in observation 3726. The first line of 
the table shows the best fit value of the carbon energy (Eq) being very low and the relative 
abundance of carbon at the preset limit value of 1000. The reduced chi-square value of the 
fit and the degree of freedom (DOF) of the fit - a measure of the signal to noise ratio of the 
fitted spectrum - indicate that the fit is reasonably good. However, since the carbon energy 
is so low and the relative abundance has to be so high to have an effect, this fit indicates 
that carbon could be neglected. In fact, attempting the fit without carbon, as shown in the 
second line of Table CD, yields a comparably good fit and does not greatly change the O and 
S initial energies, the S to O abundance ratio, or the FWHM of the lines. (We note that 
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the fits always favored line widths of around 55 eV, so after the initial fitting we have fixed 
this input variable to accelerate the calculations.) Therefore, our results indicate a negligible 
contribution of carbon to this spectrum. As we will report in a subsequent detailed report, 



we have also analyzed the XMM-Newton observations made by iBranduardi-Raymont et al. 
( 120041 l2007al ) and reach a similar conclusion. 

In Tableland Figure [2] we show the model fits to both the north and south polar auroral 
spectra (Chandra observation IDs 3726 and 4418) neglecting carbon. The inferred oxygen 
energy is near or at the largest value included in our CX atomic dataset and for sulfur it 
is within the range of about 0.7 to 1 .5 MeV/u, both in reasonable agreement with previous 
conclusions (IKharchenko et al.ll2008l ). The fits favor a sulfur to oxygen ratio greater than 
one, but we find that this is sensitive to the spectral resolution and our preliminary fits to 
XMM-Newton observations indicate smaller values of this ratio. 



4. Discussion 



If the solar wind played an important role as a source of ions leading to the CX driven 
X-ray emission near the Jovian poles, then carbon should be required to account for the 
observed spectra, because apart from the dominant components of the solar wind, hydrogen 
and helium that do not contribute to the X-ray emission, oxygen and carbon are the most 
abundant X-ray producing species (jvon Steiger et al.ll2000l ). Since our model shows no im- 
provement to the fit to Chandra observations by including carbon we conclude that ions of 
magnetospheric origin are dominant in driving this X-ray emis sion. Moreover, in other ob- 
servations, namely those detecting X-ray emission from comets (ISchwadron fc Cravensll2000l ; 
Lisse et al.ll200ll . 120051 . 120071 ; iBodewits et al.ll2007l ) clear identification of C and O lines have 
been made. In that case, the C and O (as well as N, Mg, Fe, for example) are definitely of 
solar wind origin, implying that if solar wind C contributed at Jupiter, it would likewise be 
detected. 

This conclusion also implies that solar win d carbon ions are n ot accelerated above 1 
MeV/u before precipitation into the atmosphere (ICravens et al.ll2003l ). whereas a mechanism 
for magnetospheric ions accelerates oxygen and sulfur ions to this energy or above. For 
example, as shown in Fig. [TJ the most significant flux due to carbon occurs between 300 and 
500 eV (with a peak due to the C 5+ (2p — > Is) transition at 367 eV) if given sufficiently high 
initial ion energies and abundance. However, the low carbon initial energy inferred from 
the model fit makes it unlikely that carbon ions participated in enough CX collisions in the 
high ionization stages during the precipitation to contribute appreciably. This fact is also 
reflected by the very large uncertainties accompanying the carbon abundance fits. 
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It is also interesting to note that the fits favor sulfur ion velocities (energy per mass, 
m) about half those for oxygen. This is indicative of the original sulfur ion population in 
the outer magnetosphere being mainly singly c harge d (q=l) rather than doubly charged or 



higher. That is, as suggested by I Cravens et al.l (120031 ) it is likely that a field- aligned potential 



accelerates the ions to the energies needed to produce X-rays. Then, since the total kinetic 
energy produced by a potential difference is proportional to q, the energy per mass unit 
is proportional to q/m of the original ions, suggesting a dominantly + and S + initial ion 
population with energies per mass unit for sulfur ions half those of oxygen ions. 

We also note a favoring of a greater ab undance of sulfur relati ve to oxygen in the 



present fits compared to our previous results (IKharchenko et al.ll2008l ). This comes about 
from a rather sensitive dependence of the fits on the low photon energy Chandra spectra. 
In particular the present best fit sulfur energies are generally smaller than those inferred 
previously, so that the sulfur abundance must be increased to make up for the flux at the 
lower range of X-ray energies since sulfur is the dominant source for the spectral feature 
near 300 eV. However, our preliminary fits to similar XMM-Newton observations show lower 
As for some observation dates, illustrating the sensitivity to the lower energy portion of the 
spectrum and the possible temporal variability of the ratio. Since a putative source of the 
magnetospheric sulfur and oxygen ions would be from SO x . (x— 1,2,3), further observations 
improving the reliability of the lower energy portion of the spectrum and increasing the 
coverage of it in time could provide important constraints on the relative abundance of S 
and O ions and thus their origin and the acceleration mechanism. 

Thus, additional observations would clearly be of significant help to test the CX model 
and its inputs and thereby improve our understanding of ion sources, distributions, and 
acceleration mechanism. For example, observations during other points along the solar 
activity cycle might help pick up a changing contribution from the solar wind and more and 
deeper observations would help show the range of the variability of the X-ray aurora. 
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Fig. 1. — Example synthetic spectra for the case of equal C, O, and S ion abundances and 
iVtot = 1000, Omit = 1 and £i n it = 2 MeV/u. Individual lines are marked with different 
symbols and colors for different elements (carbon: red squares; oxygen: green circles; and 
sulfur: blue triangles). Also shown (smooth lines) are convolutions of these lines with a 
Gaussian line profile of 24.5 eV FWHM and with all lines belong to a single element summed. 
The line width used in this illustration was chosen to be the same as in previous work 
(kharchenko et alJEooa ). 
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Fig. 2. — Fits to the spectra of the north and south polar X-ray auroras observed by Chandra 
(observation IDs 3726 and 4418) neglecting carbon in the CX model. The fit parameters are 
given in Table 2. 
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Table 1. Spectral fitting of the north polar spectra of Chandra observation 3726 with 
(first line) and without (second line) carbon. The columns give best fit initial ion energies 

(Ec, Eo, Es), abundances relative to oxygen (Ac, As), line width (FWHM), spectrum 
normalization (NORM), and reduced chi-square with degrees of freedom (DOF, defined as 
the number of spectral energy bins minus the number of independent inputs). All 
uncertainties are computed with A% 2 = 2.706, equivalent to 90% confidence for a single 
parameter. "*" in the errors means that the parameter is not bounded in that direction. 
The superscript outside of the bracket represents the power of 10. 



OBS ID 


E c (MeV/u) 


E (MeV/u) 


E s (MeV/u) 


A c A s 


North 3726 


0.0318:3 


2-00^.39 

1-48^.28 


0.811831 
o.8il8:S 


10001: 4.71^ 
5lf 3 


OBS ID 


FWHM (eV) 


NORM 


Red. x 2 [DOF] 




North 3726 


561L 
55 a 


(i.i^i)- 6 

(i-06l°;!)- 6 


0.84 [10] 
0.64 [13] 





a FWHM is fixed at 55 eV (see text). 



Table 2. Spectral fitting of the north and south polar spectra of Chandra observations 
3726 and 4418 without carbon. As in Table 1 all uncertainties are computed with 
Ax 2 = 2.706, equivalent to 90% confidence for a single parameter. "*" in the errors means 
that the parameter is not bounded in that direction. The superscript outside of the bracket 

represents the power of 10. 



OBS ID 


E (MeV/u) 


E s (MeV/u) 


A s 


NORM 


Red. x 2 [DOF] 


North 3726 


1-481S.28 


o.8il8:g 


°-3 


(1-0618:!)- 6 


0.64 [13] 


North 4418 


2.00l*. 37 


0.9818:1 


o c+2.6 
z -°-1.3 


(7-3l^)- 7 


1.05 [19] 


South 3726 


2-00l*. 66 


1 c i9 + * 
i - oz -0.81 


1-51^ 


(3.3l^)~ 7 


1.29 [6] 


South 4418 


2.00l*. 47 


0.6818:11 


i7l?S 


(4.2l^:8)- 7 


1.43 [11] 



